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SYSTEM AND METHOD FOR SEGMENTING TISSUE TYPES IN 
DUAL-ECHO MAGNETIC RESONANCE IMAGES 

[0001] Field of the Invention 

[0002] The present disclosure relates to image processing, and particularly to 
segmentation of dual-echo magnetic resonance images. 

BACKGROUND OF THE INVENTION 
[0003] During a magnetic resonance (MR) imaging session, the patient is placed 
inside a strong magnetic field generated by a large magnet. Magnetized protons 
within the patient, such as hydrogen atoms, align with the magnetic field produced by 
the magnet. A particular slice of the patient is exposed to radio waves that create an 
oscillating magnetic field perpendicular to a main magnetic field. The slices can be 
taken in any plane chosen by the physician or technician performing the imaging 
session. The protons in the patient's body first absorb the radio waves and then emit 
the waves by moving out of alignment with the field. As the protons return to their 
original state (before excitation), diagnostic images based upon the waves emitted by 
the patient's body are created. MR image slices can be reconstructed to provide a 3- 
dimensional picture of the region of interest. In the reconstructed image, parts of the 
body that produce a high signal are displayed as white in an MR image, while those 
with the lowest signals are displayed as black. Other body parts that have varying 
signal intensities between high and low are displayed as some shade of gray. 
[0004] Variations in image acquisition parameters affect various qualities of the 
reconstructed images generated. For example, in dual echo MR imaging, two images 
are generated for each slice, where the first and second images are generated using 
different image acquisition parameters such as different repetition times (TR) and/or 
echo times (TE). Depending on the parameters used for each image generated, one 
image may show better contrast between different certain tissues of interest, such as 
solid tissues, e.g., white and gray matter in the brain, while another image may show 
better contrast between other tissues of interest, such as fluids and solids, e.g., 
cerebrospinal fluid (CSF) and gray matter. 

[0005] Image segmentation is useful for dividing an image into meaningful regions, 
such as for distinguishing between different tissue types. Typically, each voxel of the 
image is examined and assigned a label that associates it with a region, i.e., a tissue 
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type for medical images, based on properties of the voxel, its neighbors or similarity 
to other voxels assigned to a region. The classification of voxels into tissue types 
allows for a highly intelligible display of information by the image, such as for 
educational or diagnostic purposes, for providing a guide during surgery, for defining 
boundaries of structures, for determining volumes of structures, and for monitoring 
changes in structural volumes. Multiple images taken at the same time or at different 
times may be used to study the same region of interest. Integration of information 
displayed in the multiple images is enhanced using a technique known as registration 
for correlating voxels of the multiple images. 

[0006] One method of segmentation known in the art is cluster identification in a 2- 
D scatter plot (also known as a feature space), which includes using Cartesian 
coordinates mapping the intensity of sampled voxels from a first image obtained using 
a first echo (x-axis) vs. the intensity of sample voxels from a second image obtained 
using a second echo (y-axis). Clusters of voxels in the mapping, signifying voxels 
having similar properties and likely to correspond to the same tissue type, are 
identified and assigned a color. The clusters may be formed of a tight grouping or a 
spread out diffuse grouping of mapped points. Corresponding voxels in the first and 
second image (or other registered images) are also assigned the color for a displaying 
the cluster identification segmentation within the images. Further processing may be 
performed based on the colored images, such as identification of a specific structure 
within a tissue type, and volume measurement of a structure. 
[0007] However, typically the cluster identification is performed with manual 
intervention using a qualitative analysis, which is time and resource consuming, and 
generates inconsistencies and subjectivity. 

[0008] Accordingly, there is a need for a system and method for a quantitative 
analysis of a 2-D scatter plot in which manual intervention is minimized. 

BRIEF DESCRIPTION OF THE INVENTION 
[0009] A method is provided for classifying voxels of first and second images 
generated using first and second echoes of a dual MRI scan, the first and second 
echoes corresponding to respective first and second sets of image acquisition 
parameters. The method includes the steps of providing for accessing first image data 
corresponding to the first image including an array of voxels, the first image data 
including a first associated parameter value SI for a parameter of the imaging 
associated with respective voxels of the array of voxels; providing for accessing 
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second image data corresponding to the second image including the array of voxels, 
the second image data including a second associated parameter value S2 for the 
parameter of the imaging associated with respective voxels of the array of voxels; and 
providing for determining R and 0 values associated with respective voxels of the 
array of voxels, wherein for each respective voxel R substantially equals (SI + S2 ) , 
and 0 substantially equals arctan(Sl/S2) for SI and S2 associated with the respective 
voxel. 

[0010] In another embodiment of the invention the method includes the steps of 
providing for accessing first image data corresponding to the first image including an 
array of voxels, the first image data including a first associated parameter value SI for 
a parameter of the imaging associated with respective voxels of the array of voxels; 
providing for accessing second image data corresponding to the second image 
including the array of voxels, the second image data including a second associated 
parameter value S2 for the parameter of the imaging associated with respective voxels 
of the array of voxels; providing for plotting a first plot of a first plurality of points 
having Cartesian coordinates (x,y), each point of the first plurality of points 
correlating to a respective voxel of the array of voxels, wherein for each respective 
point of the first plurality of points x is the first associated parameter value associated 
with the correlating voxel, and y is the second associated parameter value associated 
with the correlating voxel; and providing for determining R and 8 values associated 
with respective voxels of the array of voxels, wherein for each respective voxel R 
substantially equals the distance of the point correlating to the voxel as plotted in the 
first plot from an origin of the first plot, and 0 substantially equals the angle of a line 
drawn through the origin of the first plot and the point from a positive x-axis of the 
first plot. 

[0011] Steps of the method of the invention may be implemented by executing 
instructions on a processor, where the instructions are stored on a computer readable 
medium or included in a computer data signal embodied in a transmission medium. 
[0012] In another embodiment of the invention a system for performing the steps of 
the method is provided including means for accessing first image data corresponding 
to the first image including an array of voxels, the first image data including a first 
associated parameter value SI for a parameter of the imaging associated with 
respective voxels of the array of voxels; means for accessing second image data 
corresponding to the second image including the array of voxels, the second image 
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data including a second associated parameter value S2 for the parameter of the 
imaging associated with respective voxels of the array of voxels; and means for 
determining R and 8 values associated with respective voxels of the array of voxels, 
wherein for each respective voxel R substantially equals (SI 2 + S2 2 ) v \ and 0 
substantially equals arctan(Sl/S2) for SI and S2 associated with the respective voxel. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0013] FIG. 1 is a block diagram of a magnetic resonance imaging system; 
[0014] FIG. 2 is a scatter plot of image data generated, during an imaging session; 
[0015] FIG. 3 is a polar coordinate plot of the image data plotted in FIG. 2; and 
[0016] FIG. 4 is a frequency distribution histogram of the image data plotted in 
FIG. 3. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0017] Referring to FIG. 1, an exemplary magnetic resonance (MR) imaging 
system 10 is shown including an MR imager 12 and a processor assembly 14, where 
the MR imager 12 images a region of interest within a patient in an imaging session 
(also referred to as a scan) using dual-echo MR imaging, and generates image data 
including at least first and second sets of image data that correspond to imaging of the 
region of interest using first and second echoes of an echo sequence, respectively. The 
exemplary system 10 further includes a display device 16 for displaying at least one 
image that corresponds to the image data, and a user input interface 1 8 for allowing a 
user to enter information, such as data and requests, to the processing assembly 14 
and/or the imager 12. The processing assembly 14 receives and processes the 
generated image data for generating at least first and second displayable images, 
generating a 2-D scatter plot of signal intensities corresponding to the first and second 
images, generating a polar coordinate plot corresponding to the 2-D scatter plot, and 
filtering the polar coordinate plot with an edge preservation filter. 
[0018] The exemplary MR imager 12 includes a system control unit, a gradient 
amplifier system, an RF amplifier, and power supplies, and an examination chamber 
(not shown). The examining chamber includes an assembly for receiving and 
positioning a patient, a main magnet assembly including an alignment system, a 
gradient coil system and an RF transmit and receive coils. During an imaging session 
a magnetic field is applied within the examination chamber, and modified by signals 
generated by the gradient coil system and the RF transmit/receive coils. The system 
control unit controls the gradient amplifier system and the RF amplifier for applying a 
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series of radio frequency (RF) and magnetic field gradient pulses to invert or excite 
magnetic spins; and to induce, refocus, manipulate and/or encode (e.g., spatially) 
magnetic resonance. Resonance signals are received by the RF coils, demodulated, 
digitized and provided to the processing assembly 14 as 3-D image data including a 
plurality of voxels. 

[0019] An imaging sequence is determined by the RF pulses and/or the magnetic 
field gradient pulses, as controlled by the system control unit in response to operator 
entered selections. In the example provided, the image sequence selected is spin echo, 
including a series of sub-sequences including an excitation 90 degree RF pulse 
followed by first and second spaced inversion 180 RF pulses. Consecutive excitation 
pulses are spaced by repetition time (TR). Data corresponding to the first and second 
sets of echo data are acquired following the first and second inversion pulses at times 
TE1 and TE2, respectively, where the first and second sets of echo data are processed 
by the processor assembly 14 for generating first and second images, respectively, as 
described further below. 

[0020] TR, TE1 and TE2 (where the first and second inversion pulses are applied at 
TE1/2 and (TE1 + (TE2-TEl)/2 ) are operator selectable image acquisition parameters 
for acquiring images having desired characteristics, such as proton density weighted, 
T2 weighted and Tl weighted. When acquiring Tl -weighted images, short repetition 
times are used in which the Tl -weighted image is acquired as an additional sequence 
rather than a two-echo sequence. Tl -weighted images can be rapidly acquired, adding 
ittle time to another image acquisition. 

[0021] In the present example, TR, TE1 and TE2 are selected so that TR is long, 
TE1 is short and TE2 is long, for generating a first proton density weighted image in 
which tissues having a high density of hydrogen atoms are brightest, and a second T2- 
weighted image in which tissues having a long transverse relaxation time (T2), mostly 
fluids, are brightest. Typically, TR, TE1 and TE2 values for obtaining the first proton 
density weighted image and the second T2 -weighted image are substantially 
2000msec<=TR<=3000msec; 20msec<=TEK=30msec; and 

80msec<=TE2<= 100msec; where the ratio of TE2:TE1 is in the range of about 2:1 to 
5: 1 ; the ratio of TR:TE1 is in the range of about 67: 1 to 1 50: 1 ; and the ratio of 
TR:TE2 is in the range of about 20: 1 to 38: 1 . In the present example, an MR scan of a 
patient's head is performed using TR = 3000ms, TE1 = 30ms and TE2 = 80ms. 
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[0022] It is envisioned that the imaging sequence includes a series of sub- 
sequences including a first excitation pulse followed by the first inversion pulse after 
which data corresponding to the first set of data is acquired at time TE1 after the first 
inversion pulse, followed by a second excitation pulse followed by the second 
inversion pulse after which data corresponding to the second set of data is acquired at 
time TE2 after the second inversion pulse. Consecutive excitation pulses are spaced 
by repetition time TR. 

[0023] At each data acquisition data corresponding to a 2-dimensional plane (or 
slice) is acquired. Individual data acquisitions are acquired corresponding to different 
planes (or slices), so that data acquired for a series of data acquisitions corresponds to 
a 3-dimensional space. Accordingly, data corresponding to the respective first and 
second sets of image data acquired for a series of data acquisitions each correspond to 
substantially the same 3-dimensional space, where the first and second sets of image 
data are density weighted and T2 weighted, respectively. 
[0024] The processing assembly 14 includes at least one processor, such as a 
microprocessor, a CPU, a personal computer, a PDA, a hand-held computing device, 
a mainframe computer, etc. Processors of the processing assembly 14 may be in data 
communication with one another, such as by a network such as a LAN, WAN, 
intranet, internet, etc. The processing assembly 14 further includes an input port 16 
for receiving the image data. The input port 16 provides an interface between the MR 
imager 12 and the processor assembly 14. 

[0025] A variety of software modules executable by the processing assembly 14 are 
accessed by the processing assembly 14, and executed thereby for processing of the 
image data, and for determination of the condition of the region of interest. The 
software modules each include a series of programmable instructions executable on 
the processing assembly 14. The software modules may be stored on at least one 
computer readable medium (e.g., RAM, floppy, CD-ROM, flash memory, hard drive, 
etc.) or be included in a computer data signal embodied in a transmission that is 
accessible by the processing assembly 14. The at least one storage medium, and/or a 
drive associated therewith, may be external to or included within the processing 
assembly 14. The means for transmitting the signal may be partially or fully external 
to and/or included in the processing assembly 14. 

[0026] The software modules include an image generation module 30, a scatter plot 
generation module 32, a polar plot generation module 34, a histogram generation 
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module 36, and an edge preserving filter module 38. The image generation module 30 
includes an algorithm which processes the image data, including the first and second 
sets of image data, and generates a displayable 3-D image set including a plurality 
voxels Vxyz, where x, y and z indicate the a position in three dimensional space along 
the x, y and z-axes, respectively. Each voxel V xyz has associated with it first and 
second intensity signal sets SI and S2, respectively, where SI and S2 are intensity 
values that correspond to the first and second image, respectively. For spin echo 
imaging, intensity values Sl xyz and S2 xyz for a voxel V xyz are as follows, as is known 
in the art: 

(1) Sl xyz =exp(-TEl/T2 xyz )(l-exp(-TR/Tl xyz ); 

(2) S2 xyz =exp(-TE2/T2 xyz )(l-exp(-TR/Tl xyz ); 

where Tl xyz and T2 xy2 are included in the image data that corresponds to the 
voxel V xyz . 

[0027] The scatter plot generation module 32 includes an algorithm which plots for 
individual voxels V^, where the individual voxels are typically sampled voxels, on a 
Cartesian coordinate plot with S2 plotted along the X-axis and SI plotted along the Y- 
axis. An exemplary scatter plot 200 for a dual echo MRI scan is shown in FIG. 2. 
Clusters 202, 204 and 206 are shown, however a clear demarcation between the 
clusters is not apparent. Typically, a manual analysis of the scatter plot is performed 
to identify the clusters and to associate voxels in an appropriate cluster. Radial lines 
210, 212, 214 oriented at angles 01, 02 and 03 with respect to the X-axis, respectively, 
indicate estimated divisions between the clusters. The radial lines define respective 
sectors having angles 01, 02- 01, 03- 02. However, the radial line 212 is oriented so 
that a voxel 210 which belongs in cluster 206 is included with cluster 204. 
[0028] The polar plot generation module 34 includes an algorithm which calculates 
an angle 0 xyz and radius value R xyz for individual voxels V xyz , where the individual 
voxels are typically sampled voxels, in accordance with the equations: 

(3) 0xyz = arctan(Sl xy2 /S2 xy2 ) = arctan(intensity ratio) = arctan(exp((TE2- 
TE1)/T2 xyz )) 

(4) R xyz = (Sl xy2 2 + S2 xyz 2 )' /l 

Equation (3) shows that 0 correlates with T2 and indicates tissue type. The angle and 
radius values may be calculated directly from the image data received from the 



7 



29309 



imaging device 12, the 3-D displayable image data set generated by the image 
generation module 30, or the scatter plot generated by the scatter plot generation 
module 32. 

[0029] The polar plot generation module 34 may further include an algorithm that 
plots the values for 0 xy2 and R xyz on a polar coordinate plot 300, which may be provided 
to the display device 16 for display, as shown in FIG. 3. 

[0030] The histogram generation module 36 includes an algorithm for analyzing the 
frequency distribution of voxels V xyz for angles generated by the polar plot generation 
module 34 for determining the boundary or edges of clusters and the demarcation 
between adjacent clusters. The histogram generation module 36 determines the 
frequency distribution of voxels V xyz examined or sampled over values for e, i.e., the 
number N$ of voxels V*^ having a corresponding angle 0 or range of angles. Each 
cluster typically has a concentrated amount (larger frequency) of voxels, i.e., a peak, 
in its interior region relative to a diffused amount (lower frequency) of voxels, i.e., 
valleys at its boundary region. The analysis by the algorithm of the histogram 
generation module 36 recognizes the boundaries of the clusters and divisions between 
adjacent clusters by recognizing the boundary area for each cluster, which are the 
outer angles of the frequency distribution of each cluster having a minimum 
frequency for the cluster. For example, the histogram may be smoothed and the 
minima from a curve fit through the histogram would then separate the clusters. 
[0031] The histogram generation module 36 may include an algorithm for plotting 
0 versus N0, such as histogram 400 shown in FIG. 4, where 0 is plotted along the X- 
axis and N0 is plotted along the Y-axis. The width of the rectangles 402 represents a 
range of angles, and the height represents the number of voxels V xyz having an angle 
that falls within the angle range, i.e., frequency thereof. Since the angle depends on 
T2, valleys in the angle histogram indicate boundaries or edges between adjacent 
clusters, such as shown in FIG. 4 for clusters 202, 204, 206. The demarcations 
between the clusters may by determined by mathematical analysis of the 9 values for 
the voxels V xyz being examined, and/or analysis of the plotted histogram 400. In FIG. 
4, demarcations are recognized at 01b and 02b. Horizontal lines corresponding to Olb 
and 82b have been drawn on FIG. 3, showing the demarcation between the clusters 
plotted in the polar plot 300. 

[0032] A color is assigned to each cluster, so that each voxel belonging to a cluster 
is assigned the corresponding color. The image generated by the image generation 
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module 30 is redisplayed with the individual voxels displayed in the respective color 
assigned. Accordingly, the image shows segmentation in accordance with the cluster 
recognition. 

[0033] The edge preserving filter module 38 includes an algorithm which processes 
the histogram generated by the histogram generation module 34 for filtering out noise 
while preserving edges, such as by a median filter or other edge preserving noise- 
removing filter. The edge preserving filter module 38 does not change image 
parameters, such as T2 and proton density. 

[0034] The described embodiments of the present disclosure are intended to be 
illustrative rather than restrictive, and are not intended to represent every embodiment 
of the resent disclosure. Various modifications and variations can be made without 
departing from the spirit or scope of the present disclosure as set forth in the 
following claims both literally and in equivalents recognized in law. 
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